The  NBS  Daylight  Availability 
Database 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Bureau  of  Standards 
National  Engineering  Laboratory 
Center  for  Building  Technology 
Building  Physics  Division 
Gaithersburg,  Maryland  20899 


July  1984 


Prepared  for 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  California  93043 

National  Fenestration  Council 
3310  Harrison  St. 

Topeka,  KS  66611 

Naval  Facilities  Engineering  Command 
U.S.  Navy 

Washington,  DC  20390 


Directorate  of  Civil  Engineering 
U.S.  Air  Force 

j.  Washington,  DC  20330 


100 

,U56 

34-2859 

1984 


and 

Office  of  Chief  of  Engineers 
U.S.  Army 

Washington,  DC  20314 


NBSIR  84-2859 


THE  NBS  DAYLIGHT  AVAILABILITY 
DATABASE 

S.  Treado 
G Gillette 
W.  Remmert 
J.  Bean 

U S.  DEPARTMENT  OF  COMMERCE 
National  Bureau  of  Standards 
National  Engineering  Laboratory 
Center  for  Building  Technology 
Building  Physics  Division 
Gaithersburg,  Maryland  20899 

July  1984 

Prepared  for 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  California  93043 

National  Fenestration  Council 
3310  Harrison  St. 

Topeka,  KS  6661 1 

Naval  Facilities  Engineering  Command 
U.S.  Navy 

Washington,  DC  20390 

Directorate  of  Civil  Engineering 
U.S.  Air  Force 
Washington,  DC  20330 

and 

Office  of  Chief  of  Engineers 
U.S.  Army 

Washington,  DC  20314 


U.S.  DEPARTMENT  OF  COMMERCE,  Malcolm  Baldrige,  Secretary 
NATIONAL  BUREAU  OF  STANDARDS.  Ernest  Ambler.  Director 


ABSTRACT 


This  report  presents  a database  containing  hourly  measurements  of  solar 
radiation,  illumination,  sky  luminance  and  ambient  air  temperature  for  an 
entire  composite  year.  The  measurements  were  made  at  the  National  Bureau  of 
Standards,  Gaithersburg,  Maryland  (77°  west  longitude,  38°5  north  latitude). 

Both  instantaneous  hourly  and  integrated  average  hourly  measurements  are  included, 
as  are  daily,  monthly  and  annual  averages  and  totals.  For  each  measured  quantity, 
a histogram  of  the  distribution  of  the  data  is  presented  for  the  year.  The 
data  measurement,  collection  and  analysis  system  is  described. 

Key  words:  daylight  availability,  irradiance,  luminance,  sky  illuminance, 

solar  radiation,  weather  data 


iii 


TABLE  OF  CONTENTS 

Page 

ABSTRACT  iii 

LIST  OF  TABLES  v 

LIST  OF  FIGURES  v 

1.  INTRODUCTION  1 

2.  DESCRIPTION  OF  MEASUREMENTS  2 

2.1  Measured  Parameters  2 

2.2  Measurement  System  6 

2.3  Data  Analysis  and  Checking  6 

3.  DISCUSSIONS  OF  DATABASE  9 

3.1  Data  Distribution  Histograms  10 

3.2  Database  Printout  11 

4.  SUMMARY  12 

5.  REFERENCES  

13 

APPENDIX  A.  DAYLIGHT! NG  DATABASE  A-l 


iv 


LIST  OF  TABLES 


Page 

Table  1 . Sensor  Types  3 

Table  2.  List  of  Database  Parameters  4 

LIST  OF  FIGURES 

Figure  1.  NBS  Daylighting  Laboratory  14 

Figure  2.  Luminance  sensor  array  15 

Figure  3.  Diffuse  illuminance  and  irradiance  sensors  16 

Figure  4.  Vertical  illuminance  and  irradiance  sensors  17 

Figure  5.  Global  irradiance  sensor  18 

Figure  6.  Definitions  of  solar  angles  19 

Figure  7.  Data  distribution  histogram  for  zenith  luminance  (LZ)  ...  20 

Figure  8.  Data  distribution  histogram  for  north  luminance  (LN)  ....  21 

Figure  9.  Data  distribution  histogram  for  east  luminance  (LE)  22 

Figure  10.  Data  distribution  histogram  for  south  luminance  (LS)  ....  23 

Figure  11.  Data  distribution  histogram  for  west  luminance  (LW)  24 

Figure  12.  Data  distribution  histogram  for  north  irradiance  (SN)  ...  25 

Figure  13.  Data  distribution  histogram  for  south  irradiance  (SS)  ...  26 

Figure  14.  Data  distribution  histogram  for  east  irradiance  (SE)  ....  27 

Figure  15.  Data  distribution  histogram  for  west  irradiance  (SW)  ....  28 

Figure  16.  Data  distribution  histogram  for  global  illuminance  (IZ)  ..  29 

Figure  17.  Data  distribution  histogram  for  south  illuminance  (IS)  ..  30 

Figure  18.  Data  distribution  histogram  for  east  illuminance  (IE)  ...  31 

Figure  19.  Data  distribution  histogram  for  west  illuminance  (IW)  ...  32 

Figure  20.  Data  distribution  histogram  for  north  illuminance  (IN)...  33 


v 


LIST  OF  FIGURES  (Continued) 


Page 

Figure  21.  Data  distribution  histogram  for  diffuse  illuminance 

(ID)  34 

Figure  22.  Data  distribution  histogram  for  global  irradiance 

( SRT)  35 

Figure  23.  Data  distribution  histogram  for  diffuse 

irradiance  (SRD)  36 

Figure  24.  Data  distribution  histogram  for  cloud  ratio  (CR)  37 

Figure  25.  Data  distribution  histogram  for  outdoor 

temperature  (TO)  38 

Figure  26.  Data  distribution  histogram  for  direct  normal  (beam) 

illuminance  (IBN)  39 

Figure  27.  Data  distribution  histogram  for  direct  normal  (beam) 

irradiance  (SBN)  40 

Figure  28.  Data  distribution  histogram  for  atmospheric 

extinction  coefficient  El  41 

Figure  29.  Data  distribution  histogram  for  atmospheric 

extinction  coefficient  E2  42 

Figure  30.  Data  distribution  histogram  for  downward  transfer 

ratio  TR1  43 

Figure  31.  Data  distribution  histogram  for  downward  transfer 

ratio  TR2  44 


vi 


1.  INTRODUCTION 


Weather  data  has  been  collected  and  compiled  at  many  locations,  for  many 
years.  Such  information  is  used  for  a variety  of  purposes,  including  weather 
predictions,  aviation  and  agricultural  forecasts.  However,  building 
designers  need  special  types  of  weather  data  to  evaluate  the  thermal  and  energy 
performance  of  buildings.  Information  regarding  solar  radiation,  daylight  illu- 
mination and  outdoor  air  temperature  is  essential  for  the  effective  analysis  of 
building  envelopes,  the  optimum  use  of  fenestration,  and  accurate  analyses  of 
building  energy  performance.  In  some  instances,  typically  for  worst-case  or 
design  day  purposes,  single  values  are  used  for  the  pertinent  weather  parameters. 
Examples  of  this  would  be  a minimum  or  average  outdoor  air  temperature  to 
calculate  heat  transfer  through  a portion  of  a building  envelope  or  calculating 
interior  daylight  levels  assuming  an  overcast  sky.  While  single-valued  weather 
data  are  useful,  accurate  analysis  of  building  energy  performance  and  optimum 
design  of  building  components  requires  data  which  describes  the  dynamic  nature 
of  weather  parameters  over  an  annual  period.  Hourly  data  is  frequently  sufficient 
for  this  purpose. 

Most  standard  weather  stations  collect  a limited  menu  of  climatic  parameters, 
some  including  global  total  solar  irradiance  on  a horizontal  surface,  outdoor 
air  temperature,  relative  humidity,  barometric  pressure,  wind  speed  and  amount 
of  cloud  cover  [1].  These  weather  data  are  available  from  the  National  Climatic 
Center  [2]  and  various  other  sources  including  computer  time-sharing  services. 

The  standard  weather  data  sources  lack  some  detailed  information  required  for 
building  design,  such  as  data  for  radiant  and  visible  energy  incident  upon 
horizontal  and  vertical  surfaces,  luminance  distribution  of  the  sky  dome,  and 
measurements  of  cloud  conditions.  In  particular,  sky  luminance  and  illuminance 
data  are  required  for  calculating  interior  daylight  levels,  and  vertical 
irradiance  data  is  required  for  determining  fenestration  heat  transfer.  Only 
by  considerating  the  heating,  cooling  and  lighting  effects  of  fenestration 
options  can  building  designers  determine  the  most  effective  utilization  of 
fenestration. 

This  report  contains  a record  of  irradiance,  sky  illuminance,  luminance  and 
outdoor  air  temperature  measurements  over  a one-year  period.  Twenty  eight 
climate  and  solar  parameters  are  presented  on  an  hourly  basis.  Two  groups 
of  data  are  presented:  integrated  average  hourly  data  and  the  corre- 

sponding instantaneous  hourly  data.  Also  included  are  daily  totals  and 
averages,  monthly  averages  and  annual  averages.  For  each  measured  quantity,  a 
histogram  of  the  distribution  of  the  data  is  presented  for  the  year. 

The  majority  of  the  data  was  collected  during  1982,  however  entire  days  with 
missing  data  due  to  equipment  failure  have  been  replaced  with  data  from  the 
same  date  from  either  1981,  1983,  or  1984.  The  measurements  were  made  at  the  Day- 
light Laboratory,  National  Bureau  of  Standards,  Gaithersburg,  Maryland  (suburb 
of  Washington,  D.C.),  77°  west  longitude  38°5'  north  latitude  (see  figure  1.). 

The  following  sections  describe  the  measured  parameters,  data  measurement  and 
analysis  system,  followed  by  the  database. 
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2.  DESCRIPTION  OF  MEASUREMENTS 


2.1  MEASURED  PARAMETERS 

Measurements  were  made  of  diffuse  and  global  irradiance,  and  diffuse  and  global 
illuminance  on  a horizontal  surface,  global  irradiance  and  illuminance  on  vertical 
surfaces  (north-,  south-,  east-  and  west-facing),  sky  luminance  at  five  locations 
(15°  solid  angle  spots),  and  exterior  air  temperature.  Table  1 describes  the 
sensors  used  for  the  measurements.  A complete  list  of  the  measured  parameters, 
sensor  types  and  related  information  is  shown  in  table  2. 

It  should  be  noted  that  the  alphanumenic  codes  associated  with  each  measured 
parameter  do  not  correspond  to  any  particular  nomenclature  but  rather  are 
taken  directly  from  the  names  of  the  FORTRAN  variables  used  for  data  processing. 

This  provides  for  consistency  among  the  equations  in  the  text  of  the  report  and 
the  printout  of  the  database.  FORTRAN  is  limited  to  upper-case  variable  names, 
with  no  provisions  for  subscripts.  Table  2 includes  a description  of  each 
measured  parameter  and  the  associated  code. 

A general  description  of  sensor  types  and  characteristics  follows. 

a)  Sky  Luminance  - Luminance  is  the  luminous  intensity  of  a surface  in  a given 
direction  per  unit  of  projected  area  of  the  surface  as  viewed  from  that 
direction.  Luminance  was  measured  using  photovoltaic  silicon  cells, 

with  photometric  filters  and  15°  solid  angle  f ield-of-view.  One  sensor  was 
centered  on  zenith  (LZ)  and  each  of  four  other  sensors  was  centered  at  42° 
elevation,  due  north  (LN),  east  (LE),  south  (LS)  and  west  (LW).  A 
photograph  of  the  luminance  sensor  array  is  shown  in  figure  3. 

b)  Illuminance  - Illuminance  is  the  density  of  luminous  flux  (weighted  by  eye 
response)  incident  on  a surface.  Illuminance  was  measured  using  photo- 
voltaic silicon  cells,  with  photometric  filters  and  cosine  diffusers. 

Global  illuminace  (IZ)  was  measured  on  a horizontal  surface  from  an  unobstructed 
sky  hemisphere.  Diffuse  illuminance  (ID)  on  a horizontal  surface  was 
measured  using  a shading  band  (see  figure  3).  Corrections  for  the  portion 
of  sky  obstructed  by  the  shading  band  were  made  using  manufacturer's  data. 

Direct  normal  illuminance  (IBN)  was  determined  by  subtracting  diffuse  from 
global  illuminance  and  dividing  by  the  sine  of  the  solar  altitude, 


(1) 


where 


SALT  = solar  altitude 


SALT  = arc  sin  [sin  L sin  D + cos  L cos  D cos  (12-T)15] 
L = site  latitude 

D = solar  declination  = 0.4093  sin  [2JL  (J+284)] 


(3) 


(2) 
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Table  1 . Sensor  Types 


Type  A)  Luminance  - Photovoltaic  silicon  cell  with  15°  solid  angle  lens 
and  photometric  filter  (which  follows  the  standard  spectral  eye 
response  curve  [3]) 

Type  B)  Irradiance  - Photovoltaic  silicon  cell  with  cosine  corrected 
diffuser  and  radiometric  filter  (flat  response  450-950  nm 
spectral  range),  hemispherical  f ield-of-view 

Type  C)  Illuminance  - Photovolatic  silicon  cell  with  cosine  corrected 
diffuser  and  photometric  filter,  hemispherical  f ield-of-view 

Type  D)  Irradiance  - Pyranometer,  thermopile  type  (flat  response  300- 
3000  nm  spectral  range),  hemispherical  f ield-of-view 

Type  E)  Temperature  - Copper-constantan  thermocouple,  with  radiation 
shield 
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Code 

Table  2.  List  of  Database  Parameters 
Parameter  Sensor  Type* 

Measurement"*- 

Units 

Maximum 

Uncertainty 

1 

LZ 

Zenith  luminance 

A 

cdm”2 

+ 

3% 

2 

LN 

North  luminance 

A 

cdm“2 

+ 

3% 

3 

LE 

East  luminance 

A 

cdm“2 

+ 

3% 

4 

LS 

South  luminance 

A 

cdm“2 

+ 

3% 

5 

LW 

West  luminance 

A 

cdm“2 

+ 

3% 

6 

SALT 

Solar  altitude 

- 

degrees 

- 

7 

SAZ 

Solar  azimuth 

- 

degrees 

- 

8 

STIM 

Solar  time 

- 

hours 

- 

9 

SN 

North  vertical  irradiance 

B 

wm~2 

+ 

8% 

10 

ss 

South  vertical  irradiance 

B 

wm~2 

+ 

8% 

11 

SE 

East  vertical  irradiance 

B 

wm“2 

+ 

8% 

12 

SW 

West  vertical  irradiance 

B 

wm~2 

+ 

8% 

13 

IZ 

Global  horizontal  illuminance 

C 

lux 

+ 

3% 

14 

IS 

South  vertical  illuminace 

C 

lux 

+ 

3% 

15 

IE 

East  vertical  illuminance 

C 

lux 

+ 

3% 

16 

IW 

West  vertical  illuminance 

C 

lux 

+ 

3% 

17 

IN 

North  vertical  illuminance 

C 

lux 

+ 

3% 

18 

ID 

Diffuse  horizontal  illuminance 

C 

lux 

+ 

5% 

19 

SRT 

Global  horizontal  irradiance 

D 

wm”2 

+ 

4% 

20 

SRD 

Diffuse  horizontal  irradiance 

D 

wm~2 

+ 

6% 

21 

CR 

Cloud  ratio 

- 

- 

+10% 

22 

TO 

Ambient  air  temperature 

E 

°F 

+ 

PH 

O 

00 

• 

o 

23 

I BN 

Direct  normal  illuminance 

- 

lux 

+ 

8% 

24 

SBN 

Direct  normal  irradiance 

- 

wm“2 

+10% 

25 

El 

Atmospheric  extinction  coefficient  1 

- 

- 

- 

26 

E2 

Atmospheric  extinction  coefficient  2 

- 

- 

- 

27 

TR1 

Downward  transfer  ratio  1 

- 

- 

- 

28 

TR2 

Downward  transfer  ratio  2 

— 

* From 

Table 

1. 

+ cdm“2  = candelas  per  square  meter 
Wm“2  = watts  per  square  meter 
lux  = lumens  per  square  meter 
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J = Julian  Date  1 _<  J 365 
T = Solar  time. 

Illuminances  on  vertical  planes  (north,  east,  south  and  west)  were  measured 
using  an  array  shielded  from  ground  reflections  (see  figure  4). 

c)  Irradiance  - Irradiance  is  the  density  of  radiant  flux  (equally  weighted  for 

all  wavelengths)  incident  on  a surface.  Global  irradiance  (SRT)  on  a horizontal 
surface  was  measured  using  a thermal-type  pyranometer  exposed  to  the  unobstructed 
sky  dome  (see  figure  5).  Diffuse  irradiance  (SRD)  was  measured  with  a similar 
Instrument  and  a shading  band.  Corrections  for  the  portion  of  the  sky 
obstructed  by  the  shading  band  were  made  using  manufacturer's  data.  Direct 
normal  irradiance  (SBN)  was  also  determined  by  subtracting  diffuse  from 
global  irradiance  and  dividing  by  the  sine  of  the  solar  altitude. 


SBN  = 


SRT  - SRD 
sin( SALT) 


(4) 


d)  Air  temperature  - Outdoor  air  temperature  was  measured  using  a type-T 
(copper-constantan)  thermocouple,  with  a radiation  shield. 

The  database  also  contains  the  derived  quantities  of  solar  altitude,  solar 
azimuth,  solar  time,  atmospheric  extinction  coefficient  and  downward 
transfer  ratio.  Figure  6 shows  the  definitions  of  the  solar  angles. 

e)  Solar  altitude  - Calculated  using  equation  2. 

f)  Solar  azimuth  - Calculated  using  equation  5a,  b,  or  c. 


W = ( 12-T)  x 15  (5) 

If  cos  W < tan  D/tan  L 

SAZ  = arcsin  (cos  D sin  W/cos  SALT)  (5a) 

If  cos  W > tan  D/tan  L 

SAZ  =>  it  - arcsin  (cos  D sin  W/cos  SALT)  (5b) 

If  SAZ  < 0,  SAZ  = SAZ  + 2tt  (5c  > 

Solar  azimuth  is  zero  due  north,  90°  due  east,  180°  due  south,  and  270°  due 
west . 

g)  Solar  time  - Calculated  from  equation  6. 

STIM  = TS  + (1  if  daylight  time)  + ET  - 2/15  (6) 
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(7) 


where  TS  = standard  time  (decimal  hours) 

ET  = equation  of  time 

= 0.170  sin  ( J-80)  -0.129  sin  2 (j-8) 

373  355 

h)  Cloud  ratio  = The  ratio  of  diffuse  to  global  irradiance. 

CR  = SRD  (8) 

SRT 

i)  Atmospheric  extinction  coefficent  - Two  versions  (E1,E2)  were  determined  from 
measurements,  based  on  direct  illuminance  or  irradiance,  as  follows: 

El  (direct  illuminance) 

El  = - sin  (SALT)  Jin  IBN  (9) 

IXT 

where  IXT  = 127,500  { 1 + 0.034  cos  -2JL  (j-2)  } (10) 


E2  (direct  irradiance) 

E2  = - sin  (SALT)  Jin  J3BN 

SXT 

where  SXT  = 1353  I 1 + 0.034  cos  2JL  (J-2)  } 

365 


(11) 

(12) 


j)  Downward  transfer  ratio  - Two  versions  (TR1  and  TR2)  were  computed,  one 
from  irradiance  and  one  from  illuminance.  The  downward  transfer  ratio  is 
defined  as  the  portion  of  solar  energy  scattered  out  of  the  direct  beam 
reaching  the  ground  [7].  Therefore  one  minus  the  downward  transfer  ratio 
is  reflected  back  into  space,  or  absorbed  by  the  atmosphere. 

TR1  (illuminance) 

TR1  = ID  (13) 

IXT  sin( SALT) [l-EXP(-El/sin( SALT))] 

TR2  (irradiance) 

TR2  = _ SRD,  (14) 

SXT  sin( SALT) [l-EXP(-E2/sin( SALT))] 

2.2  MEASUREMENT  SYSTEM 


The  data  acquisition  system  consisted  of  a microcomputer-controlled  data  logger 
which  monitored  the  measurement  sensors.  All  sensors  were  configured  to  provide 
voltage  inputs  to  the  data  logger,  which  performed  analog-to-digital  conversion 
and  voltage  measurement.  The  data  logger  also  provided  temperature  compensation 
for  the  thermocouple  readings.  Under  control  of  the  microcomputer,  sensor 
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readings  were  continuously  and  sequentially  scanned.  At  the  beginning  of 
each  hour,  a data  file  for  that  hour  was  assembled  consisting  of  the  current 
reading  for  each  sensor  (instantaneous  value)  and  the  average  of  all  the  readings 
for  each  sensor  over  the  previous  one  hour  period  (average  value).  Each  average 
value  is  calculated  from  fifteen  readings  during  the  hour.  The  letter  "A" 
preceding  a measured  parameter  name  indicates  an  average  value. 

Data  were  stored  on  0.133  m (5.25  in.)  floppy  disks,  and  transferred  to  a larger 
computer  for  subsequent  processing,  analysis  and  storage. 


2.3  DATA  ANALYSIS  AND  CHECKING 

The  database  contains  measurements  of  eighteen  parameters,  and  ten  derived 
quantities.  Over  the  period  of  a year,  this  amounts  to  approximately  2.5 
million  readings,  including  average  and  instantaneous  values.  Each  annual  aver- 
age parameter  is  a result  of  131,400  readings,  and  each  instantaneous  parameter 
is  represented  by  8760  readings.  The  initial  raw  data  files  were  processed  to 
check  for  errors  and  inconsistencies.  A careful  review  of  the  data  was  conducted 
to  identify  erroneous  data,  however,  data  were  only  removed  from  the  database 
if  a clear  justification,  as  described  below  could  be  made  for  such  removal. 
Consequently,  if  a datum  passed  all  the  checking  routines  and  was  not  associated 
with  any  known  measurement  problems,  it  was  included  in  the  database.  Even 
though  the  checking  routines  were  considered  a successful  technique  for  controlling 
erroneous  data,  from  time  to  time  questionable  data  still  appeared  and  were 
flagged  accordingly.  The  shear  volume  of  data  processed  could  have  resulted 
in  a small  portion  of  inaccurate  data  remaining  in  the  database  undetected. 

Potential  inaccuracies  can  be  grouped  into  two  general  categories.  The  first 
category  includes  what  could  be  described  as  measurement  errors.  In  this 
category  would  fall  calibration  errors,  sensor  non-linearity  and  errors  in 
voltage  measurement  or  analog-to-digital  conversion.  The  second  category 
includes  measurements  which  are  accurate,  but  not  measuring  the  parameter  of 
interest  properly.  Prime  examples  of  this  would  be  a sky  luminance  measurement 
when  the  solar  disk  is  visible  in  the  field  of  view  of  the  luminance  sensor,  or 
when  snow  is  present  on  sensors.  In  these  cases,  the  measurement  system  may  be 
working  properly,  but  the  reading  would  not  be  representative  of  sky  luminance 
as  assumed. 

The  first  stage  of  data  checking  included  adjusting  each  reading  for  any  zero 
offset,  setting  to  zero  any  slight  non-zero  readings  known  to  be  zero  (i.e., 
irradiance  after  dark),  and  setting  diffuse  irradiance  and  illuminance  equal  to 
total  irradiance  and  illuminance,  respectively,  if  the  respective  diffuse  com- 
ponent reading  exceeded  the  reading  for  total.  This  second  correction  is  occa- 
sionally necessary  on  overcast  days  due  to  the  average  nature  of  the  shading 
band  correction  factor,  but  is  usually  a small  effect. 

The  second  stage  of  data  checking  involved  comparing  each  reading  to  an 
acceptable  range  for  each  parameter,  and  intercomparing  related  readings  as  a 
group.  Maximum  and  minimum  levels  were  established  for  each  parameter,  based 
on  the  expected  range  of  theoretically  possible  readings.  Readings  outside  of 
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this  range  were  flagged  with  an  asterisk  (*)  and  set  to  999999.  For  the  lumi- 
nance channels,  a reading  in  excess  of  15,000  cdm~2  Was  assumed  to  be  due  to 
interference  from  the  solar  disk,  and  therefore  flagged  with  an  asterisk  but 
not  set  to  999999  unless  the  reading  exceeded  75,000  cdm~2 . Three  groups  of 
parameters  were  assembled  for  intercomparison.  The  five  luminance  sensors 
were  one  group,  the  four  vertical  irradiance  sensors  the  second  group  and  the 
four  vertical  illuminance  sensors  the  third  group.  Within  each  group  an  average 
of  each  hour’s  readings  was  computed,  and  if  any  individual  reading  differed 
from  the  average  by  more  than  a factor  of  ten  (or  1/10)  that  reading  was  flagged 
with  a plus  sign  (+)  and  set  to  999999. 

The  averaged  hourly  data  should  be  more  accurate  than  the  instantaneous  hourly 
data  for  two  reasons.  First,  any  random  errors  (electronic  noise,  analog  to 
digital  conversion,  etc.)  should  cancel  out  when  multiple  readings  are  averaged. 
Second,  an  occasional  false  reading  due  to  intermittent  malfunctions  or  condi- 
tions will  have  a small  effect  on  the  average  reading.  Both  instantaneous  and 
average  readings  are  of  value.  In  many  daylighting  applications,  the  instantan- 
eous conditions  are  of  interest.  However,  a one  hour  average  is  nearly  "instan- 
taneous" in  relation  to  an  annual  file  containing  8760  hourly  readings.  Since 
many  energy  calculations  are  done  on  an  hourly  basis,  perhaps  hourly  averaged 
daylighting  data  are  appropriate  to  use  for  these  applications. 

When  data  were  unavailable  for  1982,  due  to  equipment  down  time,  data  were 
substituted  from  1983,  1984  or  1981  for  the  same  julian  date.  While  this  changes 
the  distribution  of  the  data  and  affects  the  average  slightly,  the  usefulness 
of  the  annual  data  is  not  compromised,  since  1982  should  not  be  considered  a 
special  year.  Thus,  the  database  reflects  reasonable  conditions  for  a possible 
year,  but  not  necessarily  a typical  year. 
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3.  DISCUSSION  OF  DATABASE 


The  primary  objective  of  the  publication  of  this  database  is  to  provide  a 
comprehensive  set  of  reference,  concurrent  measurements  of  illuminance, 
irradiance  and  sky  luminance  for  a one-year  period.  These  data  are  critically 
needed  for  building  design  purposes  and  daylighting  analyses,  but  are  virtually 
unavailable . 

All  or  portions  of  the  database  might  be  used  as  driving  parameters  for  energy 
or  daylighting  calculations,  or  for  other  purposes  such  as  sky  and  solar  model- 
ing. Detailed  analysis  of  the  relationships  between  the  measured  parameters, 
or  their  dynamic  behavior,  may  enable  new  correlations  to  be  developed  for 
improved  sky  and  solar  modeling.  Since  the  annual  database  contains  such  a 
vast  quantity  of  information,  it  would  be  useful  to  analyze  the  hourly  data  and 
determine  means  to  characterize  the  measurements  in  a more  general  manner. 

Two  methods  were  used  to  characterize  and  summarize  the  measurements.  One 
method  involved  computing  averages  for  each  of  the  parameters  on  a daily,  a 
monthly  and  an  annual  basis.  The  second  method  consisted  of  determining  the 
mean  and  standard  deviation  of  each  parameter  and  plotting  a histogram  showing 
the  distribution  of  the  data  over  the  year. 

The  computed  averages  of  each  parameter  provide  a means  for  making  an  assessment 
of  the  general  magnitude  of  that  parameter,  and  allow  comparison  of  different 
day's  data.  For  example,  a very  bright  day  will  have  a much  higher  average 
illuminance  than  a dark  overcast  day.  In  addition  to  the  daily  average  for 
each  parameter,  which  is  computed  using  only  the  non-zero  readings,  the  daily 
total  of  each  measured  parameter  was  tabulated.  This  quantity  indicates  the 
total  of  all  readings  for  each  parameter,  thereby  including  the  effects  of 
longer  or  shorter  daylight  periods  throughout  the  year.  Note  that  averages 
based  on  24  hours  can  easily  be  computed  by  dividing  each  daily  total  by  24. 

Monthly  averages  were  also  calculated  for  each  parameter.  The  monthly  average  is 
defined  as  the  average  daily  total  for  the  month  (not  the  average  hourly 
reading).  Annual  averages  were  calculated  for  each  parameter,  being  equal  to 
the  average  daily  total  for  the  year  and  the  hourly  average  of  non-zero  read- 
ings. Exceptions  to  the  above  criteria  were  required  for  the  atmospheric 
extinction  coefficients  and  downward  transfer  ratios,  since  a total  atmospheric 
extinction  coefficient  for  the  day  would  have  no  meaning.  Here  only  hourly 
averages  for  each  month  and  the  annual  average  were  tabulated. 

While  the  daily,  monthly  and  annual  averages  provide  much  useful  information, 
they  do  not  give  an  indication  of  the  variability  or  distribution  of  the  values 
of  each  parameter.  An  effective  way  of  presenting  the  distribution  of  each 
parameter  is  through  a histogram  of  the  readings  showing  the  percentage  of  the 
entire  annual  data  set  falling  into  various  intervals.  Such  histograms  were 
developed  from  the  data  to  show  the  spread  of  the  readings  for  each  parameter, 
along  with  the  probability  of  an  individual  reading  falling  within  a certain 
range.  The  mean  and  standard  deviation  of  the  annual  set  of  data  for  each 
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parameter  were  also  calculated.  For  this  usage,  the  standard  deviations  do  not 
represent  an  uncertainty,  but  rather  indicate  the  variability  of  each  parameter 
throughout  the  year. 

Before  discussing  the  individual  histograms,  a general  examination  of  the 
dynamic  nature  of  the  database  parameters  is  appropriate.  The  behavior  of  each 
of  the  parameters  through  the  year  is  governed  by  the  superposition  of  several 
cycles.  The  irradiance,  illuminance  and  luminance  data  follow  a strong  diurnal 
(daily)  cycle  in  which  maximum  values  typically  occur  near  a certain  set  of 
solar  angles.  For  example,  energy  incident  on  east-facing  vertical  surfaces  is 
typically  maximum  in  the  morning,  while  horizontal  and  south-facing  surfaces 
typically  receive  maximum  incident  energy  at  solar  noon.  However,  this  does 
not  mean  that  south  vertical  irradiance  is  always  maximum  at  solar  noon,  since 
considerable  hourly  variations  can  occur  due  to  a randomly  fluctuating  cycle 
superimposed  on  the  diurnal  cycle.  This  random  cycle  is  a function  of  the 
local  weather  conditions  as  influenced  by  cloud  cover,  atmospheric  water  vapor 
and  other  factors.  A third  significant  cycle  is  due  to  seasonal  effects.  Many 
geographical  locations  exhibit  observable  differences  in  available  solar  radia- 
tion or  cloud  conditions  on  a seasonal  basis,  whereby  more  overcast  days  occur 
in  winter  than  in  summer  (or  vice  versa).  Outdoor  air  temperature  follows  a 
strong  seasonal  cycle  and  a strong  diurnal  cycle.  Direct  normal  irradiance  and 
illuminance  are  less  sensitive  to  a daily  cycle,  except  for  the  effects  of  low 
solar  elevations. 

3.1  DATA  DISTRIBUTION  HISTOGRAMS 


A pair  of  histograms  is  presented  for  each  parameter,  the  first  based  on 
instantaneous  readings,  the  second  based  on  average  readings.  Figures  7 through 
11  present  the  distribution  of  readings  for  each  of  the  sky  luminance  sensors. 
Sky  luminance  shows  a wide  range  in  the  data,  but  a fairly  narrow  band  of 
probable  intensity.  Readings  greater  than  15,000  cdm~2 t which  are  due  to 
the  solar  disk  being  present  in  the  luminance  measurement  area,  are  excluded 
from  these  figures. 

Figures  12  through  15  present  the  distribution  of  vertical  irradiance  readings. 
As  would  be  expected,  north  vertical  irradiance  is  significantly  lower  than 
the  other  vertical  planes  averaging  only  one  third  of  south  vertical  irradiance. 
East  and  west  vertical  irradiances  are  about  two  thirds  of  the  south  vertical 
readings . 

Figures  16  through  21  present  the  distribution  of  the  illuminance  readings. 
Considerable  variation  is  apparent  for  all  surfaces,  graphically  illustrating 
the  wide  dynamic  range  of  daylight  availability.  Diffuse  horizontal  illumin- 
ance and  east,  south  and  west  vertical  illuminance  range  from  66  to  50  percent 
of  global  horizontal  illuminance,  while  north  vertical  is  less  than  25  percent. 

The  distributions  global  horizontal  irradiance  and  diffuse  horizontal 
irradiance  are  presented  in  figures  22  and  23,  respectively.  Diffuse  irradi- 
ance average  about  50  percent  of  global irradiance . The  ratio  of  diffuse  to 
global irradiance , cloud  ratio,  is  shown  in  figure  24.  The  average  cloud  ratio 
was  68  percent,  but  a wide  range  occurred. 
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Figure  25  shows  the  distribution  of  outdoor  air  temperatures.  The  mean  air 
temperature  was  10.6  C (51.1°F). 

Direct  normal  illuminance  and  irradiance  are  presented  in  figures  26  and  27, 
respectively.  They  exhibit  a fairly  even  distribution.  The  computed  atmospheric 
extinction  coefficients  are  shown  in  figures  28  and  29.  Over  half  of  the  read- 
ings were  less  than  one  or  greater  than  5.75.  The  downward  transfer  ratios  are 
presented  in  figures  30  and  31.  TR1  (illuminance)  averages  36  percent,  while 
TR2  (irradiance)  averages  31  percent. 

3.2  DATABASE  PRINTOUT 


Appendix  A contains  the  printout  of  the  daylighting  database.  Each  parameter 
is  labeled  with  an  alpha  code  as  listed  in  table  2.  One  day's  worth  of  data  is 
on  each  page,  with  adjacent  pages  containing  hourly  average  data  for  the  day, 
and  instantaneous  hourly  data  for  the  day.  The  julian  date  and  calendar  date 
are  listed  at  the  top  of  each  page.  Each  page  is  divided  into  two  sections 
with  twelve  hour's  data  in  each  section.  Hours  are  indicated  by  column  headings. 
The  units  of  measurement  are  listed,  along  with  daily  totals  and  non-zero  aver- 
ages for  each  parameters.  The  last  page  of  the  appendix  contains  the  daily  or 
hourly  averages  (as  appropriate)  on  a monthly  and  annual  basis,  for  each 
parameter . 

Missing  or  defective  data  were  represented  by  999999.  Uncertain  or  questionable 
data  were  flagged  with  an  asterisk  (*)  or  plus  sign  (+) . No  attempt  was  made 
to  synthesize  replacement  data  for  missing  or  bad  data,  because  the  intent  of 
this  paper  is  to  provide  measured  data,  rather  than  calculated.  However,  the 
user  can  use  whatever  means  they  feel  is  appropriate  to  fill  in  for  missing 
data.  Procedures  such  as  interpolating  between  available  data  points  [1],  or 
calculating  based  on  other  parameters  [4]  or  models  [5,6]  would  be  useful  for 
that  purpose. 
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4 . SUMMARY 


An  annual  database  is  presented  including  irradiance,  illuminance,  sky  luminance, 
sky  condition  and  temperature  measurements.  The  database  is  intended  for  use 
in  daylighting  design  and  calculations,  building  energy  analyses,  and  sky  and 
solar  modeling.  The  measurement  system  is  described.  The  data  is  analyzed  and 
processed  including  examination  of  reading  averages,  distribution  and  variability 
over  a one  year  period. 
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Figure  1.  NBS  Daylighting  Laboratory 
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Figure  2.  Luminance  sensor  array 
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Figure  3.  Diffuse  illuminance  and  irradiance  sensors 
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Figure  4. 


Vertical  illuminance  and  irradiance  sensors 
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Figure  5.  Global  irradiance  sensor 
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Data  distribution  histogram  for  zenith  luminance  (LZ) 
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Data  distribution  histogram  for  east  luminance  (LE) 
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Figure  10.  Data  distribution  histogram  for  south  luminance  (LS) 
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Figure  11.  Data  distribution  histogram  for  west  luminance  (LW) 
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Figure  12.  Data  distribution  histogram  for  north  irradiance  (SN) 
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Figure  13.  Daita  distribution  histogram  for  south  irradiance  (SS) 
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Figure  14.  Data  distribution  histogram  for  east  irradiance  (SE) 
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Figure  15.  Data  distribution  histogram  for  west  irradiance  (SW) 
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Figure  17.  Data  distribution  histogram  for  south  illuminance  (IS) 
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Figure  18.  Data  distribution  histogram  for  east  illuminance  (IE) 
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Figure  19.  Data  distribution  histogram  for  west  illuminance  (IW) 
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Figure  20.  Data  distribution  histogram  for  north  illuminance  (IN) 
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Figure  24.  Data  distribution  histogram  for  cloud  ratio  (CR) 
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Figure  25.  Data  distribution  histogram  for  outdoor 
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Figure  26.  Data  distribution  histogram  for  direct  normal  (beam) 
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Figure  27.  Data  distribution  histogram  for  direct  normal  (beam) 
irradiance  (SBN) 


40 


AV.  =1 

EACH  ♦ lEPRfc  SENT  S 35  CM  SERVATISP.S 


II33LE  3F 

NUMBER 

OF 

NTERVA. 

OBSERVATIONS 

0.0 

330 

********** 

0.5 

1002 

• ***************************>  ****************  * 

1.0 

532 

**************** 

1.5 

245 

******* 

2.0 

125 

**** 

2.5 

76 

*** 

3.0 

46 

** 

3.5 

26 

* 

4.0 

11 

* 

4.5 

4 

* 

5.0 

5 

* 

5.5 

2 

* 

5.0 

1056 

*********************** »»»»;»« * 

N 

Aif  • El 

4064 

mean 

2.153 

ME  DI  AN 

0.820 

THE  AN 

2.052 

ST  jE  V 

2.365 

SEfliAN 

0.037 

MAX 

5.000 

min 

0.010 

33 

6.000 

01 

0.400 

IN.  El 

EACH  * REPRc  SEKIS  35  33 SER VA T I 3 RS 


MIDDLE  3F 

NUMBER 

OF 

INTERVA. 

OBSERVATIONS 

0.0 

303 

********* 

0.5 

1585 

************* 

1.0 

431 

************* 

1.5 

225 

******* 

2.0 

133 

**** 

2.5 

91 

*** 

3.0 

61 

** 

3.5 

29 

* 

4.0 

21 

* 

4.5 

10 

* 

5.0 

13 

• 

5.5 

7 

• 

6.0 

1143 

************* 

IN.  El 

H 

4052 

MEAN 

2.305 

MEDIAN 

0.850 

TMcAN 

2.220 

STD:  V 

2.428 

SEME  AN 

0.036 

MAX 

5.000 

MIN 

0.010 

03 

5.000 

01 

0.390 

Figure  28.  Data  distribution  histogram  for  atmospheric 
extinction  coefficient  El 


41 


A/.  E2 


£ A CH  9 RE 

PRt  SEMS 

33  (JlSERXATIuNS 

MIDDLE  DF 

NUMBER 

DF 

I NIT ERVAL 

LBSE  RV  AT  T 343 

0.0 

323 

99999999999 

0.5 

1 497 

99999999999999999999999999999  9999999 

1.0 

53  e 

•99999999999999999 

1.5 

2 41 

999999999 

2.0 

124 

#9999 

2.5 

73 

999 

3.0 

E7 

999 

3.5 

37 

99 

4.0 

35 

99 

4.5 

13 

9 

5.0 

1C 

9 

5.5 

10 

9 

6.0 

11C7 

99999  999  99  99999999999999999999999999. 

N 

AX.  E2 

4075 

MEA  i 

2.277 

MEOI AM 

0.900 

TNEA  N 

2.189 

STDEV 

2.39b 

SEHSAN 

3.036 

MAX 

6.000 

MIN 

3.010 

03 

5.000 

01 

0.400 

in.  e 2 

EACH  * REPRESENTS  35  □JSERVATI  uNS 


IID3LE  DF 

NUMBER 

OF 

MTERXA_ 

OBSERVATION* 

0.0 

2 77 

99999999 

0.5 

1542 

999999999  99999999999999999999999999999 9 9*  9999 

1.0 

411 

*99999999999 

1.5 

149 

99999 

2.0 

115 

9999 

2.5 

92 

999 

3.0 

66 

99 

3.5 

45 

99 

4.0 

32 

9 

4.5 

35 

9 

5.0 

21 

9 

5.5 

11 

9 

6.0 

1287 

999999999999999999999999999999999999* 

N 

IN.  E2 
4083 

MEAN 

2.523 

MEDIAN 

0.930 

THEAM 

2.462 

STDEV 

2.504 

SEME  AN 

0.039 

MAX 

6.000 

MIN 

0.010 

03 

5.000 

01 

3.390 

Figure  29.  Data  distribution  histogram  for  atmospheric 
extinction  coefficient  E2 
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Figure  30. 
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Figure  31.  Data  distribution  histogram  for  downward  transfer 
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